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Abstract 

Radiation by solid or fluid bodies can be characterized by resonance modes. 
They are complex, as well as resonance frequencies, because of the energy loss 
due to radiation. For ducts, they can be computed from the knowledge of the 
radiation impedance matrix. For the case of a flanged duct of finite length ra- 
diating on one side in an infinite medium, the expression of this matrix was 
given by Zorumski, using a decomposition in duct modes. In order to calculate 
the resonance frequencies, the formulation used in Zorumski's theory must be 
modified as it is not valid for complex frequencies. The analytical development 
of the Green's function in free space used by Zorumski depends on the integrals 
of Bessel functions which become divergent for complex frequencies. This pa- 
per proposes first a development of the Green's function which is valid for all 
frequencies. Results are applied to the calculation of the complex resonance fre- 
quencies of a flanged duct, by using a formulation of the internal pressure based 
upon cascade impedance matrices. Several series of resonance modes are found, 
each series being shown to be related to a dominant duct mode. Influence of 
higher order duct modes and the results for several fluid densities is presented 
and discussed. 

Keywords: Acoustics, radiation impedance, cylindrical pipe, resonance 
frequencies 
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1. Introduction 

For the problem of duct radiation, many works have been done concerning 
the calculation of radiation for a given duct mode, but to the author's knowkedge 
no study has been done concerning the resonance modes. We will treat the 
problem, in order to get a better insight of the coupling of the duct and the 
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surrounding space by radiation. We choose the case of a duct with infinite 
flange, because of its relative simplicity. One difficulty is due to the fact that 
resonance frequencies (and modes) are complex, because radiation is a form 
of dissipation. Notice that in the literature, resonance modes are also called 
eigenmodes: they must be distinguished from the duct modes used in the present 
paper for the purpose of the calculation. 

Green's functions are widely used in many physical situations and notably 
in acoustics, for the calculation of the pressure field radiated by physical sources 
{e.g. speakers, musical instruments, vibrating structures,...). For instance, the 
solution given by Rayleigh [ij to the classical problem of a plane piston radiating 
into an infinite flange involves the Green's function in free space. Zorumski [2| 
extended this result to know the radiation of a semi-infinite fianged duct in the 
form of matrix impedance, giving the coupling between duct modes and used the 
Sonine's infinite integral ( Ref. [3'|, p. 416, Eq.4) to develop the Green's function 
in free space. In Rcfs [4, 5], formulations based on the Zorumski's method of 
this radiation impedance are obtained for a larger class of problems. However, 
the development fails when the frequency becomes complex: the corresponding 
infinite integral, involving a Bessel function whose argument is a product of the 
frequency and the dummy argument, becomes divergent for complex frequen- 
cies. In many studies, the Zorumski's radiation matrix is used as a boundary 
condition at the end of the duct in order to calculate input impedances, length 
corrections or reflection coeflicients (see e.g. Refs. or Q). It is worth noting 
that complex resonance frequencies can occur in various situations {e.g. dissi- 
pative fluid, radiation, complex impedance wall boundary conditions such as in 
Refs. [1] or [§|). In section [2] of this paper, we present a new expression for the 
Green's function in free space for complex frequencies. In section [U an appli- 
cation of this result is devoted to the determination of the complex resonance 
frequencies of a cylindrical duct, closed at its input, considering the influence 
of higher order duct modes. In the same secion, some results are given and dis- 
cussed. For this purpose, the internal Green's function is previously calculated 
in section [3] with a method of cascade impedance. 



2. Calculations of Green's function for the Helmholtz equation in free 
space 

Many studies on sound radiation by cylindrical ducts can be found in the 
literature. For the case of an infinite flange, Norris and Sheng [l3| or Nomura 
[Hj used a Green's function integral to find an appropriate formulation for the 
external field. We can also cite the classical work by Levine and Schwinger 



12| for the case of an unhanged pipe. Zorumski [2] extended the results for the 
planar mode to obtain a multimodal radiation impedance which is a combination 
of the duct modes present in the duct. In this section, these calculations are 
briefly recalled, exhibiting the diflBculty related to complex frequencies. Thus, a 
new analytical formula for the Green's function, valid for a dissipative problem, 
is presented. 
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2.1. Zorumski's radiation impedance 

We consider the radiation of sound into an infinite half space from a circular 
duct (with radius b and length L), with an infinite fiange at zq = (the index 
corresponds to the cross section Sq at the end of the duct) and we have chosen 
to work with circular coordinates where the vector r is denoted (z, r, 9) as shown 
by Fig. m 




z 



«o=0 



Figure 1: Schema and coordinates of the duct. 

The acoustic pressure in the infinite medium (z > 0) is given by a Helmholtz 
integral (the time factor exp(— icji) is omitted throughout this paper): 

1*2,77 i*h ikh 

P(r) = -^y^ roviro,9o) — drodeo, (1) 

where h = [r^ + Tq — 2rrQCOs{6 — 6*0) + z^]^, p the ambient density and the 
wavenumber k = uj/c (with oj the circular frequency and c the speed of sound). 
The pressure p and the velocity v inside the duct (z < 0) are expressed as a 
series of duct eigen modes, so in z = 0: 

p(r, e,z^O)=pc^Y.T. ^mn(kr)e'"''Pmn, (2) 

7n n 

v{r, e,z^0)^cJ2Yl V-mJfcOe^'Knn, (3) 

m n 

where iprnn{kr)e^^''^ is the transverse function for the mode mn with ipmn{kr) = 
Jo(kr)/Nmn- The Xmn are the eigenvalues, solutions of J'(Xmnkb) = 0, and the 
norm Nmn is chosen similarly to that used by Zorumski [2|. 
Substituting Eq. ^ into Eq. (P) gives the pressure for z > in terms of the 
modal velocity amplitudes Vmn- 

i()nr f'^'^ e^^'^ 

p{r,e,z) = ~-^Y.T.^^-^ / ro— ^„„(A:ro)drod0o. (4) 

Zorumski expressed the free space Green's function in equation ^ in terms 
of a Sonine's infinite integral (Q, p. 416, Eq. 4) and wrote for z in the 
expression of h: 

ikh j-oo 

——^k t{t^ -l)-2j^{Tkh)dT. (5) 

Jo 
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Next, he introduced a concept of " generalized radiation impedance matrix Zray" 
for a semi-infinite duct with an infinite flange to describe the relation between 
the modal pressure and velocity amplitudes: 

oo 

Pjnji ^ ^ ^innl^mh (^) 
1=1 

where m, I and n are, respectively, the orders of circumferential, radial incident 
and reflected modes. The element Zmni of the radiation impedance matrix 
gives the contribution of the velocity mode ml to the pressure mode mn (for 
reasons of symmetry, the coupling is possible only for a duct mode with the same 
azimuthal dependance) . The expression of the radiation impedance is obtained 
as: 

Zmnl = -i t{t^ -iy^D,n,i{T,k)D„rl{T,k)dT, (7) 

Jo 

with, for a hard wall condition: 

D^,.ir,k)^kb^-^t^^^^^ (8) 



mn 



2.2. Green's function for the Helmholtz equation in free space for complex fre- 
quencies 



The following asymptotic form (see Ref. [13|, Eq. 9.2.1, p. 364) occurs when 
V is fixed and \n\ — >■ oo: 

Mn) = J^[cos{k - --k) + el^(«)IO(|«|-i)], (9) 

V TTK I 4 

with |argK| < tt (in this paper, the real part and imaginary part are repre- 
sented, respectively, by the symbols 3? and 3). As a consequence, for r — > oo 
with fc € C we have Jo{Tkh) — >■ oo when 3(fc) ^ 0, thus relation ([5]) and the 
radiation impedance ([7]) given by Zorumski are divergent integrals for all non 
real frequencies. 

In order to have a Green's function for the Helmholtz equation in free space 
valid for complex frequencies, we use another form of the Sonine's infinite in- 
tegral to develop this Green's function, expressed by Watson [3] (p. 416, Eq. 
4): 



^ikh 



t{t^ ~k^y2jQ{Th)dT. (10) 



This integral remains convergent even for k complex. A difficulty occurs since k 
is a branch point of the square root. For a time factor exp(— jwt), the integration 
path on the real axis must remain below k. However the complex resonance 
frequencies u = ck have a negative imaginary part (see explanation in section 
2]), so the previous formula must be adapted because of the branch cut. The 
integration path below k is classically deformed, as shown in Fig. [2] 
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Figure 2: Deformation of the integration contour. 

Now, the integral in Eq. (flUl) is written as: 

ikh fMk)\ roo 

— = Jo{Th)T{T^-kY^dT+h{k,h)+Ic+l2{k,h)+ Jo{Th)T{T^-e)—2dT, 

" Jo "'|K(A:)| 

(11) 

with 

hik, h) ^ hik, h)^- JoiTh)T{T^ - k^y^dr 

J|!R(fe)| 

and Ic is zero according to Jordan's lemma. 



After calculations (similar to those developed by Morse and Feshbach [14 1 
p. 410), the following five cases can be distinguished: 
i) 3?(/c) < and 5(fc) > 



— = MrhMr' -er-dT. (12) 



ii) 3fi(fc) < and 3(fc) < 

-7— = -i I M'rh)^===dT + Jo{Th)—===dT 

" Jo ^/k'^-T^ Jisft(fc)l \/T--k'^ 

l-k 

-2 / Jo(r/i) ^^^ rfr. (13) 

J|!R(fe)| VT^ - k- 

iii) 3fi(fc) > and 3(fc) > 

ikh /•|K(fc)| /-oo 

— = J,{rh)^==dT + J,(rh)^==dT. (14) 

iv) 5R(fc) > and 3(fc) = 



Jo{kh)kV2e{l - i) + i / jQ{Th)—===dT 

\Jk^ — 



Q 



Jo(T/i)^===dr, (15) 
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with e << 1. 

v) 3?(fc) > and < 



gjfc/i /■|3?(fe)| 



= +i Mrh) =dT+ / Mrh) ^ dT 





-2 r Jo(r/^)^=L=dT. (16) 

Contrary to the original Zorumski's formulation, these results involve con- 
vergent integrals when the frequency is complex. 

2.3. New formulation of generalized impedance of a flanged circular duct for 
real frequencies 

Initially, the previous result will be checked for the real case, in order to 
compare the original Zorumski's result, noted in Eq. ([7]), and that obtained 
using the expansion of exp{ikh) / h for the case (iv) where = 0. Eq. (|15p 
leads to the following results: 



T 



-ik[D„in{k)D,nl{k)k\/2e{l-i)+i I -j===Drnn{T)D„il{T)dT 



where 



/ ^=L==A„n(r)D„,i(r)dr], (17) 



UmnKT) = — — ^ . (18) 



The radiation impedance for the planar mode (m = n = with I — 0) 
with Zorumski's formulation ([7]) and formulation ([T7|) (with e = 10~^) are very 
similarly and thus, confirms the validity of formula (jl7p with the identical com- 
putational cost. This formula is used in Ref. (l5| to calculate an approximation 
of the reflection coefficient and of the length correction, taking into account the 
effect of the higher order duct modes below the first cut-off frequency. 
The comparison between complex Zorumski's formulation for the planar mode 
(m=n=0 and 1=0) and Rayleigh's radiation impedance of a flanged plane piston 
confirm the validity of Zqoo for all the frequencies (see subsection 14. 1|) , because 
Rayleigh's radiation impedance is by definition the same quantity as Zqqq. 
In what follows, we show the interest of the radiation impedance valid for com- 
plex frequencies when calculating the complex resonance frequencies of a flanged 
finite length duct terminating in a Zorumski's radiation condition. In a first in- 
stance, the internal Green's function must be determined, with a method of 
cascade impedances. 
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3. Calculation of the finite flanged duct internal Green's function 
with a method of cascade impedances 



3.1. Definition oj the internal Green's Junction 

We search for the internal Green's function G{M, M',uj) at a point M (r, 6, z) 
with a source in M'{r' , 9' , z'), satisfying: 

(Am + fc2)G(M, M\iu) = --^(5(r - r')S{e - e')S{z ~ z'), (19) 

Zirr 

with dnG(M, M' ,uj) = on the walls. It is classically (see e.g. Morse and 
Feshbach lj|) expanded in a series of duct modes (the boundary conditions for 
the variable z will be given later on): 

oo oo 

(r' ,0')gmn{z,z',uj), (20) 

rn—O n— 

where 0mn('", ^) = V'mn(^)e*™^ is the transverse function for the mode mn and 
gmn{z, z',uj) is the longitudinal function for the mode mn. In what follows, the 
dependance in w of (7m„(z, z', oj) is omitted for simplicity. It is worth noting 
that now we have the transverse function with respect to r and not kr, as per 
Zorumski [2]. Therefore, ipmnir) = '^"'^ where the Xmn are solutions of 
J'rai^ninh) = 0, with A„i„ = ^nin/b. The Tmn are the (n + 1)*'* zeros of the first 
derivative of the Bessel's function Jm- The following norm N„m is chosen (Ref. 
with respect to h and not kh): 



Nmn = by/nJ 1 - ■j^—^Jm{>^Tnnb). (21) 

The transverse modes 4)mn{f,0) satisfy: 

(Ai + A^j0„„(r,0) = O, 

with Aj_ = IMr-^ + thus: Aj_0™„(r,6') = -A^„0™(r, 6*). Studies 

such as Refs. [6| or [7| show that a small finite number of terms is necessary for a 
numerical estimate of the summation and it can be truncated: Mm is defined as 
the number of circumferential modes m and Nn as the number of radial modes 
n. Introducing the previous expression ()20p in ()19p gives: 

E E(7rT + klJ^mnir, e)4>mn{r\ e')gmn{z, z') = -—5{r - r')5{e - 9')d{z - z'), (22) 

m— n—0 

where fc^„ = fc^ _ A^„. Then, multiplying the left and right sides of ([22)) by 
(/<^,;(r, ^^) and integrating the resulting equality on the surface S, the orthogo- 
nality condition (/^ <pmn4'm'i^^ ~ Srnm'Sni), Icads to: 

{-g^ + kmn)9nin{z,z') ^ -S{z ~ z'). (23) 
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With a conservative Neuraann or Dirichlet boundary condition at the ex- 
tremities of the duct, the resonance wavenumbers can be easily computed. But 
with a dissipative boundary condition hke those occurring for the sound radi- 
ation of a flanged cylinder, there is no simple analytical solution. Therefore, 
in the next section, the elements gmniz, z') are calculated with a method of 
cascade impedances presented in Refs. [6|,[l6| or |17|. 



3.2. Presentation of the method of cascade impedances 

The calculation of the duct impedance at an abscissa zi with respect to 
another abscissa Z2 is based on the following transfer matrix relationship (see 
e.g. Ref. 11^): 

cosh(ifc„„(2;2 - zi)) Zc,,nn sinh(iA:,„„(2;2 - ^^i)) \ 
Mr ^ I sinh(zfc^„(z2-zi)) eosh(zfc_(z2 - z,)) ' ^^^^^ ^'^'^ 

. ... ' 

is an element of the diagonal matrix of characteristic impedance Zc: 

_ ^PH 

Moreover, matrices formulation is chosen: $(r, 6) is a column vector constituted 
by the M^Nn elements (f>mn, verifying Jg^{r,9)^^ {r,9)dS — I, P{z) is a 
column vector constituted by the MmNn elements Pmn and 'V{z) is a column 
vector constituted by the MmNn elements Vmn- Thus, relation ([M)) is now 
written: 

• P(zi) ZeS \ / P(Z2) 

V(zi) i I Z-^S C M V(Z2) 



(25) 



where C and S are diagonal matrices constituted by the elements cosh(ifcm„(z2 — 
Z\)) and s\Tih{ikmn{z2 — zi)), respectively. 

The transfer matrix formulation (1251) is now transformed into an impedance ma- 



trix formulation. The calculation, presented in Ref. 17[, is recalled in Appendix 
A. This gives the following matrix equation: 

V P(Z2) y I Z21 -Z22 ) \ V(Z2) ) ' ^ ' 

where Zn = ZcS-^C, Z12 = ZcS-^, Z21 = ZcS-^ and Z22 = ZcS-^C. 



3.3. Calculation of P(z') 

First, the pressure vector at the source position z' is calculated. For this 
purpose, we calculate a right-side matrix impedance Z+(2;') at z' with respect 
to Zray, a Icft-sidc matrix impedance Z (z') at z' with respect to Ze and fi- 
nally the connection between these two matrices at z' , the abscissa of the source. 
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step 1: Right-side matrix impedance T/^iTl) 
Let us denote Cg and Sg the diagonal matrix constituted by the elements 
cosh(ifcmTi^2) and sinh(ifcmn^2), respectively, where Z2 (see Fig. [T]) is the distance 
between the abscissa z' and the extremity {z — 0). The radiation impedance 
Zray, Constituted by the elements Zrani calculated in the previous section, ver- 
ifies P(z' + h) = ZrayV(z' + ^a)- Thus Eq. (US]) leads to: 

V(z' + h) = (Zray + T^llT^'LlX^ {z') , (27) 

and, with Eqs. ^ and p7|. to: 

P(z') = ZiiV(z') - Zi2(Zray + Z22 ) " ' Z21 V(z') • 

With Z+ verifying P(z') ^ T^{z')N{z'), we have finally: 

Z+(z') = Zii - Zi2(Zray + Z22)"'Z21, (28) 

or: 

Z+(z') = ZeS-^Cg - ZcS-l[Z-lZray + K^C^Y^K^. (29) 



Step 2: Left-side matrix impedance Z~(z') 

We choose a Neumann condition for z = —L (thus V(z' — l\) = 0). Here, Ce 
and Se are the diagonal matrices constituted by the elements cosh(iA;mTt/i) and 
sinh(ifcm„Zi), respectively, li (see Fig. [1]) is the distance between the point 
z = —L and the point z' . The impedance Ze is calculated at z = —L. 
Relation ([^ is written for the present case as: 

P(z') = ZeS, iV(z' - ^i) - ZcS-iCeV(z')- (30) 

Thus, with the Neumann condition in z = ~L and with Eq. (j30p . using P{z') = 
Z~(z')V(z'), the following result is obtained: 

Z-(z') = -ZcS-iCe. (31) 

Step 3: Connection between the impedance matrices Z+ and Z^ 
at z' 

Let us denote P^^i^') = [gmn{z, z')]z=z'±f Using the continuity of the Green's 
function at z = z', leads to when e — >■ 0: 

Pmn{z ) — Pmn{z ) = Pmn{z ), (32) 

Integrating relation ((23)) on an interval of width 2e between z' + £ and z' — e 
gives: 

rz'+e q2 

/ (^ + kl^n)9mn{z, z')dz = -1 (33) 

and, with the pressure continuity, when e — >■ 0: 

5.P+„(z')-5.P_(z')=-l, (34) 
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with 5zP^„(z') = [dzg,nn{z,z')]^=^,±^. 

Euler's dimensionless equation j;^Vmn{z') ~ ^^dzPmn{z') implies dzPmn{z') = 
ikVmniz'), thus, using Eq. (IMl) : 

V+Jz') - V-Jz') = -1. (35) 

We introduce a column vector W of MmNn lines whose elements equal 1. Equa- 
tion pSI) may be expressed as: 

iZ+{z'))-'P+iz') - {Z-iz'))-'p-{z') = -Iw, 

IK 

and using Eq. (|32]) . the pressure at the source is written as follows: 

P(z') ^ -1 [(Z+(z'))-' - (Z-(z'))-']"' W. (36) 

3.4- Expression of the function g(z,z') 

We introduce a column vector g(z, z') constituted by the MmNn elements 
gmn{z,z'). They are two possible configurations with respect to the relative 
positions of receiver at z and source at z' : 

First configuration: z > z' 

Let Ir = z — z' he the distance between the receiver and the source (see Fig. [Ij, 
Ci, the diagonal matrix constituted by the elements cosh.{ikmnlr) , and Si^ the 
diagonal matrix constituted by the elements sinh(ifcmn/r)- Relation ([24|) gives 
for z > z': 

g(z,z') = Ci^P(z')-ZcS,,V(z'), 
then, with P(z') = Z+(z')Y{z'): 

g(z, z') = Si^ [S,;iC,^ - Zc(Z+(z'))-'] V{z'). (37) 

Second configuration: z' > z 

Let li = z' — z he the distance between the receiver and the source(see Fig. [T]), 
Cij the diagonal matrix constituted by the elements cosh{ikmnh), and Sij the 
diagonal matrix constituted by the elements sinh(iA:,„„/;). Relation p4)) gives 
for z' > z: 

g(z,z') = Ci,P(0') + ZcSi.V(z'), 
then, with P(z') = Z-(z')V(z'): 

g(z, z') - Si, [Si;iC,, + Ze(Z-(z'))-'] P{z'). (38) 

Finally, the Green's function of a finite duct with an infinite flange is given as: 

G{M, M') = *(r, e)^{r' ,e')g{z, z'). (39) 
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4. Application to complex resonance frequencies of a flanged, flnite 
length duct 

Resonances of a flanged, finite length duct are interesting as they contain 
important information about the coupling between internal and external fluids. 
Their calculation is based on the fact that the internal pressure becomes inflnite 
at each resonance. Newton's method is used to compute the zeros of the inverse 
of the pressure. Since the resonances of a dissipative problem are complex, a 
complex formulation of the impedance radiation is needed. As a time depen- 
dence exp(— zwi) has been chosen, the imaginary part needs to be negative for 
resonance frequencies in order to ensure that the amplitude remains bounded 
for all times t. Using the integrals and for 5R(/c) > and 3(fc) < 0, 
Zmni becomes: 

Zmnl = -ik[i / —===DmniT)DmliT)dT+ / , „ ,„ -Pmn('r)A 

Jo VK - T"' "'|SR(fe)l V - 

\ I 9 ,„ D7nn{T)Dml{T)dT]. 

This expression is used as a radiation condition Zray at the end of the finite 
length duct. 



4.1. Resonance wavenumbers for the planar mode without influence of higher 
order duct modes 

In a first instance, we consider the resonance wavenumbers of the planar 
duct mode (m=n=0) without the influence of higher order duct modes (1=0). 
With radiation, the j*'' resonance wavenumbers of duct mode mn are denoted 
fc^„ J, and denoted fc^j„ without radiation. 

In order to validate the complex formulation of the radiation impedance, we 
compare the resonances obtained with radiation impedance given by relation 
(j40p for m = n = and 1 = 0, with the resonances obtained with the radiation 
impedance of a flanged plane piston given by Rayleigh's formulation as (see Ref. 
(3 P-1458): 

Zo pc[l - ^Ji{2kb) ^ ^Si{2kb)] (41) 

The radiation impedance for m=n=l=0 calculated with relation (j40l) and 
that of a flanged plane piston (j4ip are identical. Thus, resonance wavenumbers 
calculated with these two formulations are so identical, as observed in Fig. [3] 

It is worth noting in Table [T] that without radiation, the resonance frequen- 
cies are those obtained by the usual longitudinal resonances of a cylindrical duct 
with one side "closed" and the other side "open" (Neumann/Dirichlet problem) 
for: 

_ {2j + l)n 
'^00 — or I J — ■•• 
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Figure 3: Evolution of resonance wavenumbers fcpQ ^ depending on wavenumber k (with b = 
0.1m and L = Im) with the radiation impedance defined by Eq. I|40|l for m = n = and 
I = (+) and with the radiation impedance of a flanged plane piston (Eq. I4ip (□). 



We can observe in Table [T] of Appendix B that the real part of resonance fre- 
quency decreases when radiation is taken into account: this is normal behavior 
because the reactive effect of radiation can roughly be described as an increase 
in the duct length. The semi-infinite duct length correction is estimated by the 
following formula (see Ref. [loj with the modification given in private commu- 
nication): 

AT" 1-1- 

^ = 0-82159 ' . (42) 

1.2949 ^ V 1.2949 ' 

The difference between the real part of the j first resonance wavenumbers /cqq ^ of 
a finite radiating cylindrical duct (below the first cut off wavenumber fccutoi = 
3.83/6) and those estimated using the length correction defined by Eq. (|42p 
(fegp = (2j -I- l)7r/[2(L -I- Ai)]), tends to zero when the ratio L/b increases. 
The values agree well (< 1% of difference) for L/b> 3 with several values of L. 

4-2. Resonance wavenumbers with influence of higher order duct modes 

The principal interest of the complex Zorumski's formulation is that we 
can observe the influence of higher order duct modes (also denoted by H.M 
afterward). In this paper, we taken into account only the axisymmetric modes 
(m = 0) (but the non-axisymmetric are very easy to calculate with the same 
method). Figure S] shows resonance wavenumbers /cqq ^, /cq^ fcpa r- We observe 
three series of resonances. Each series starts at the cutoff frequency of a duct 
mode. The first one corresponds to a domination of the planar duct mode, the 
second one to a domination of duct mode 01 and the third one to a domination 
of duct mode 02 (see Appendix B, the Green's function profile is plotted around 
/cqq and fcoi)- ^'^S- [S]shows the infiuence of the two first higher order duct modes 
{m = 0,1 = 1 and m = 0, / = 2) on the resonance wavenumbers of the first series. 
It is worth noting that below the first cut off wavenumber kcutm = 3.83/6, only 
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one higher order duct mode is sufficient to accurately describe the resonances 
(see some values in Table [T]); between the first and second cut off wavenumber 
kcut02 = 7.02/5, only two higher order duct modes are enough and similarly to 
the higher order duct modes: between the n"* and (n+ 1)*'' cut off, only {n+ 1) 
higher order duct modes are enough. 



00 AA . 



Figure 4: Resonance wavenumbers Aiqq ^ (+), fcgj^ ^ (a), fcg2 ^ (*) for L = Im and b = 0.1m. 
Figure [B. 81 in Appendix C shows the mode profiles around the two wavenumbers indicated by 



an arrow. 



10 20 30 



40 50 60 70 
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Figure 5: Resonance wavenumbers /Cqq ^ of the first series {m = 0,n = 0) with an influence of 
(1 = 0: *), 1 (1 = 1: A) and 2 (l = 2:'+) H.M for L = Im and b = 0.1m 
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J 




fc^o^ with H.M fcgp^ with 1 H.M 




1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 



1.5708 
4.712 
7.854 
10.996 
14.137 
17.279 
20.42 
23.562 
26.7 
29.85 
32.99 
36.13 
39.27 



1.449 - 0.0095i 1.451 - 0.0096i 

4.369-0.077i 4.375-0.0776i 

7.333-0. 182i 7.345-0.183i 

10.336-0. 296i 10.345-0.299i 
13.365-0.412i 13.4-0.415i 
16.409-0. 526i 16.45-0. 53i 

19.463-0.645i 19.523-0.644i 

22.523-0. 773i 22.608-0.76i 

25.59 -0.924i 25.707-0.881i 

28.674-1. 12i 28.824-l.Oli 

31.822-1.408i 31.965-1. 152i 

35.377-1. 711i 35.149-l.317i 

38.747-1.409i 38.395-1.579i 



Table 1: Values of the j first resonance wavenumbers without radiation (k^g) and with radi- 
ation (^QQ J,) for and 1 H.M, with b = 0.1m and L = Im. 

4-3. Evolution of the j first resonance wavenumbers with respect to radiation 

In the present section, we show the evolution of the j first resonance wavenum- 
bers when only the planar duct mode propagates with respect to radiation and 
as shown in the previous section, we take into account the effect of one higher 
order duct mode. For this purpose, we introduce a multiplicative coefficient on 
the radiation impedance. Physically, this coefficient rjp can be regarded as the 
ratio between the external fluid density pext and the internal fluid density pmt, 
such as: 

Vp = , (43) 

Pint 

the density of the external fluid pext varying from a vacuum to water density, 
the sound celerity, 340m. s~^, remaining constant. 

Figures[6]and[7]show that when the parameter rjp increases, a Neumann/Neumann 
problem is obtained: the real part tends to jtt/L and the imaginary part tends 
to zero. This behavior corresponds to a system without losses, the external fluid 
becoming a perfectly reflective surface. 

Fig. [7] shows that the absolute value of the imaginary part of the resonance 
frequency goes through a maximum, corresponding to the maximum of the ra- 
diation. Similarly results have been observed for non planar modes. So, we can 
conclude that energy radiation losses evolves with the densities of internal and 
external fluid and goes through a maximum for a specific densities ratio. 



5. Conclusion 

A development of the Green's function for the Helmholtz equation in a free 
space valid for complex frequencies is possible and leads to a new formula in Zo- 



14 



Figure 6: Evolution of the real part of the j first longitudinal resonance wavenumbers fep, 
with respect to rip, for L = Im and b = 0.1m. 




-0.7' ' ' ^ 

10"^ 10"^ 10° 10 



Figure 7: Evolution of the imaginary part of the j first longitudinal resonance wavenumbers 
fcpQ ^ with respect to rjp, for L = Im and b = 0.1m. 

rumski's radiation impedance. The interest has been shown for an appUcation 
example, dedicated to the calculations of the complex resonance frequencies of 
a radiating flanged cylindrical duct. It has been shown that length correction 
calculated for a semi-infinite duct is a good estimate of a finite duct radiation 
when the ratio L/b > 2 and for frequencies sufficiently below the first cut off 
frequency. The study of the influence of higher order duct modes has shown 
that below the first cut off frequency, only one higher order duct mode is needed 
to accurately describe the influence of the external fluid on the resonances and 
between the n*'' and (n + 1)*'' cut off, only {n + 1) higher order duct modes are 
enough. In the last part, it has been observed a maximum of radiation for a 
specific densities ratio. In the future, it will be interesting to use a BEM method 
to study more complicated geometries and to observe the resonances with an 
experimental method. This work is a first step to study the relation between 
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radiation and several parameters in order to optimize geometry for minimizing 
(e.g. for noise pollution) or maximizing the sound radiation (e.g. for wind in- 
struments). 
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Appendix A. Matricial calculation 

The main steps required to obtain the results of section 13.21 are presented 
here (see Ref. (T^). 

For a cylinder, the general solutions at point zi can be described with respect 
to the values of P and V at point Z2 such as described by relation ([M)) . For all 
the modes mn, the following matrix problem is obtained: 

P(zi) =CP(z2) + ZeSV(z2), (A.l) 

V(zi) = Z-iSP(z2) + CV(z2), (A.2) 

C being a diagonal matrix constituted by the elements cosh(ifc,„„(z2 ~ -^^i)) et 
S a diagonal matrix constituted by the elements sinh(ifcm„(z2 — zi)). Eq. (jA.2p 
implies: 

P(z2) - ZeS-iV(zi) - ZcS-iCV(z2). (A.3) 

Introducing Eq. (jA.3|) in Eq. ()A.ip and using the commutativity of the 
diagonal matrices, we obtain: 

P(zi) = ZeS-lCV(zi) - [ZeS-^CC - ZeS]V(z2), 

with 

ZcS-^CC-ZcS = ZcS-^(I + SS) - ZcS 
= ZqS ZqS — ZjjS 
— ZcS ^, 

thus: 

P(zi) = ZeS-iCV(zi) - ZcS-iV(z2). (A.4) 
Therefore, with Eqs. (|A.3p and (|A.4p . we obtain the relation (f^. 
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Appendix B. Green's function profile in the duct around resonance 
frequencies 

Figure IB. 81 shows that around the resonance frequency /cqq, the profile of 
Green's function corresponds to the profile of the planar duct mode even if duct 
mode 01 is propagating and similarly around the resonance frequency kQ^, the 
profile of Green's function corresponds to the profile of the first non planar 
duct mode even if planar duct mode is propagating. The same comportment is 
observed around other resonance frequencies. Therefore, it is worth noting that 
each series observed in Fig. |4] corresponds to a predominant duct mode even 
if other duct modes are propagating. Notice that the evanescent duct modes 
exist mainly near to the source (at z = —0.5): this clearly appears on the upper 
figure. 



m=0 and n=3,l=3 with k=44.63-2.246i 




m=0 and n=3,l=3 with k=44.65-0.235i 




z 



Figure B.8: Profile in the duct of the real part of Green's function around k^^ with 3 H.M 
(upper figure) and around k^-^ with 3 H.M (lower figure) 
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